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Abstract. A thermodynamics for systems at a stationary state is formulated. It is 
based upon the assumption of the existence of local equilibrium in phase space which 
enables one to interpret the probability density and its conjugated nonequilibrium 
chemical potential as mesoscopic thermodynamic variables. The probability current 
is obtained from the entropy production related to the probability diffusion process 
and leads to the formulation of the Fokkcr-Planck equation. For the case of a gas of 
Brownian particles under steady flow in the dilute and concentrated regimes we derive 
nonequilibrium equations of state. 



PACS numbers: 05.70.Ln, 82.70.Dd, 83.60.Fg 



Mesoscopic thermodynamics of stationary states 



2 



1. Introduction 

Systems driven outside equilibrium by the intervention of an external driving force 
frequently show a peculiar behaviour related to their relevant quantities not observed 
when they are close to equilibrium. The dependence of thermodynamic quantities 
on external gradients |I]-|SI and the long-range |3] and aging |H] behaviours of the 
correlation functions are examples illustrating this very common characteristic 
The explanation of those peculiarities has been the subject of intense activity in recent 
years. A common feature observed is that a description performed by simply proceeding 
through the extension of equilibrium concepts to an out-of-equilibrium situation does 
not always lead to a correct characterization of the state of the system |12[ I13j. Such 
an approach can only be justified when the system is close to equilibrium in a local 
equilibrium state in which the relevant variables are the locally-conserved hydrodynamic 
fields ^3]. In any other circumstance taking place when other variables become relevant, 
one has to propose a more general scenario able to provide a complete description of 
the system. 

Our purpose in this paper is precisely to analyze one case presenting those typical 
features: a 'gas' of Brownian particles under a shear flow fH]- We will propose a new 
thermodynamic theory aimed at the characterization of the system at a stationary state. 
The theory is based upon a broader interpretation of the concept of local equilibrium 
at the mesoscale |in]-|IH] which enables us to analyze irreversible processes taking place 
in time scales where fluctuations become manifest. We will derive the nonequilibrium 
equation of state for the pressure tensor as a function of the shear rate in the dilute 
and concentrated regimes finding the non-analytic dependences observed in simple and 
complex fluids P El E] • 

The paper is organized as follows. In Sec. 2, we derive the Fokker-Planck equation 
describing the dynamics of a dilute gas of Brownian particles beyond the point particle 
approximation and obtain the hydrodynamic equations from the hierarchy of evolution 
equations for the moments of the probability distribution. In section 3, expressions 
of the mobility and transport coefficients depending on frequency and shear rate are 
derived from hydrodynamics through the Faxen theorem. Non-equilibrium equations of 
state valid in the dilute limit are proposed in Section 4 whereas Section 5 is devoted to 
the case of higher concentration. Finally, in the discussion section we summarize our 
main results. 

2. Thermodynamics of a 'gas' of Brovi^nian particles in external flow 

We consider a 'gas' of Brownian particles in contact with a heat bath and subjected 
to conditions creating a velocity field Vq = VQ{f,t). The dynamic description of the 
particles is accomplished by means of the single-particle probability density P{f,u,t), 
where r is the position of a particle and u its instantaneous velocity. This probability 
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is normalized to the number of particles and obeys the conservation law 

|p + V.(5P) = -|..(P..). (1) 

where Pv^i represents an unknown probability current in -u-space and V = 

The explicit expression of the current Pv^, can be found by analyzing the dissipation 
taking place in w-space, or equivalently, the entropy production. To this end, we will 
propose a mesoscopic thermodynamic description based on the assumption of local 
equilibrium in phase space leading us to formulate the Gibbs equation 

Tpds = p6e — / fiSPdu, (2) 



in which the entropy depends on the probability distribution function instead of on the 
density of mass. Here p{r,t), s{f,t), e{f,t) and p{r,t) are the density, the entropy, the 
internal energy and the nonequilibrium chemical potential per unit of mass and T the 
temperature, assumed constant. The density field is defined as 

p{f,t) = my P{u,f,t)du, (3) 

where m is the mass of the particle, and the velocity field v{f, t) is given through the 
first moment 



pv{f,t) = 1^ J uP{u,f,t)du. (4) 
The Gibbs equation formulated in Q must be compatible with the Gibbs entropy 



f P 

s{t) = -ks / caln—du + sie. (5) 

J -He 

P 



Here fc^ is the Boltzmann constant, = — the mass fraction of the Brownian particles, 
sie the local equilibrium entropy and Pie the probabihty density of the local equilibrium 
reference state given by ^H] 

P,^ = e^''^''=-^("-"^°^'l, (6) 

where pie is the chemical potential at local equilibrium, and ^{u — VqY the kinetic energy 
per unit of mass corresponding to the relative motion of the particles with respect to 
the flow. Making variations with respect to P in Eq. (jS)) and comparing the result with 
(jSJ we infer the value of the chemical potential 

p=^\nP+hu-vor. (7) 
m 2 

The entropy production follows by taking the time derivative of Eq. Q and by 
using Eqs. (jT}, Q, (jHI), and the corresponding equation for [T^. After integrating by 
parts and assuming that the probability current vanishes at the boundaries, one obtains 



Pv^ ■ -^du -- j J -Vy-iu- vof\ - ^ y • Fdu, (8) 



where -F = is the force per unit of mass exerted by the flow on the particle. In 
deriving Eq. ©, we have also used the balance equation for the internal energy e. 
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assumed isothermal conditions, neglected the energy dissipation arising from viscous 
heating [T3] and defined ^ = ■^+v-V. The first contribution to the entropy production 
(jHl) is related to the diffusion process in -u-space. The second originates from the relative 
current J = {u — v)P, whereas the third, involving the current Jq = {u — vq)P comes 
from the dependence on time of the external velocity field vo{r,t) 

According to the Onsager theory ^3], we may now establish couplings between 
fiuxes and forces of the same tensorial order appearing in Eq. In particular, the 
diffusion current in u is given by 



Pv,, 



ou 



[u - vo) -Vvo + C- FP, 



(9) 



where we have introduced the coefficients ^, e and ( which may, in general, be 
functions of position and time |T^. Their dependence on position accounts for spatial 
inhomogeneities whereas the dependence on time introduces memory effects CHI 1201 • 
By substituting now Eq. (jU)) into (P), one arrives at the Fokker-Planck equation 



d_ 
di 



P + V • (uP) 



d_ 

du 



m ou 



d 






C- FP 


du 



,(10) 



where we have introduced the friction tensor [3 related to ^ through (3 = ^ + e- Vvq. 

In this equation which describes the dynamics of the system beyond the point 
particle approximation, the inertial and surface forces the bath exerts on the particle 

— * 

during its motion are associated to the terms Q ■ F and e • Vvq, respectively. 

From the Fokker-Planck equation, we can directly derive the evolution equations 
for the moments of the distribution which correspond to the hydrodynamical fields 
In particular, the mass conservation equation 

l^-V-M (11) 

can be obtained by taking the time derivative of equation (jS)), using (fTUI) and integrating 

by parts. By following a similar procedure, from Eq. Q we may derive the momentum 

conservation law 

dv ~-* ~-* 

+ =-p(3-{v-Vo)+pC-F, 



dt 



->k 



in which we have introduced the kinetic part of the pressure tensor P (r, t) 

P {r,t)=m {u — v){u — v)P{u, f, t)du, 
whose evolution equation is 

2kBT ^ 



m 



-Pi ■ 



(12) 



(13) 



(14) 



Here, the exponent s refers to the symmetric part and the matrix of relaxation times r 
is defined by 



1 ^1-1 
(5 + Vv + -{V -v)! 



(15) 
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In Eq. (|14j) . we have neglected the contribution of the third centered moment and 
higher order moments related to the Burnett and super-Burnett approximations since 
their characteristic relaxation times are much smaller than those appearing in Eqs. (jllj) - 
()14|) 13'. It is important to stress that Eq. ()14|) constitutes a general expression from 

which one may obtain the constitutive equation for P which may incorporate memory 
effects. 

At times t ^ the particle enters the diffusion regime governed by the 

corresponding Smoluchowski equation for p{r,t) jT3]. In this regime, Eq. ()14|) yields 

P ~ -^pl - [fe + r/^) • ViTol', (16) 

where we have assumed that for sufficiently long times, v ~ + 0(V Inp) and have 
been defined the Brownian ff^ and hydrodynamic ffj^ viscosities 

rfo = — — pP and ifj, = — — p/3 ■ e. (17) 
m m 

The former originates from the stresses caused by the Brownian particle in its motion 

whereas the latter is a consequence of the finite size of the particle which introduces 

inertial effects. Eq. (fTBj) will constitute the starting point for the derivation of the 

nonequilibrium equations of state of the gas. 

3. The generalized Faxen theorem 

The explicit form of the coefficients /3, e and can be obtained by making use of the 
Faxen theorem [21] giving the force exerted by the fluid on a particle of finite size, which 
has been generalized to the case of a shear flow [22] • For simplicty, we will neglect the 

contribution due to Brownian motion in Eqs. (fT^ and (fTB)) . and assume that {3 = f31, 

e = el and (f = CI- 

By substituting Eq. into p2j] after neglecting Brownian contributions, 

simplifying for p and taking the Fourier transform of the resulting equation, one obtains 



v{u;) = G{u) 



k T 

l3vo(uj) H pP~^eV'^Vo{uj) - iujQvo^uj) 

m 



where we have defined the propagator of the velocity G = {—iu + /5)~^, and v{lo) = 
(27r)~i / ve~^'^^dt. The quantities ^0(1^) and V'^vq{uj) appearing in Eq. (fT5|l can be 
expressed as a linear combination of the averages of the unperturbed velocity over the 
surface and the volume of the particle defined respectively as 

^^(u;) = / vo{^)dS = voico) + ^V^voi^), (19) 



Ana^ J ' ' ' ' 3! 

and 

^^(^) = 1^ / vo{^)df= vo{u) + ^V%iuj), (20) 
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where a is the radius of the particle and the derivatives of the velocity are evaluated at 
the center of mass. Using these relations, Eq. (fTHjl can be written as 



ViUJ 



m a-^ 
"^15 , 



m 



3^ 3^. 

2/9+ 

' -B + -Ciuo 



Vo [iO), 



(21) 



which coincides with Faxen's theorem obtained from hydrodynamics fIT] if we make the 
following identifications 

1 



1 + aa^ + -a'^al 
9 



1 

e = -- 



m 



QksT 



59 



1 + 2aa^ + -^(^^(^l I ) 



(22) 
(23) 



and C = where f3o is the Stokes friction coefficient, a^j = \J is the inverse of the 
penetration length of the perturbation, v the kinematic viscosity of the solvent and Pp 
and pf are the particle and fluid densities. 

The case of a stationary flow Vq = Vo{r) and zero- frequency can also be analyzed by 
simply considering the average over the surface of the particle since in this case inertial 

effects vanish. The mobility tensor f3 can also in this case be compared with the one 
calculated from hydrodynamics |221 



One also obtains 



1 m 



-a 



5 



(24) 



(25) 



where = 



^ is the inverse of the penetration length of the perturbation related to 



the shear flow, with 7 the shear rate and M a matrix whose elements depend on the 
boundary conditions on the surface of the particle |23j . 



4. Nonequilibrium equations of state 

Once the expresions of the transport coefficients have been found, we will proceed to 
derive the equations of state for the different elements of the pressure tensor as functions 
of the shear rate 7 in the zero-frequency case. 

Since the gas of Brownian particles constitutes a fluid with a viscosity coefficient 
proportional to the diffusiveness Dq = with r]f being the viscosity of the heat 

bath, in this case we will introduce the inverse penetration length = where we 

have replaced the kinematic viscosity p of the heat bath by the diffusivity Dq given by 
the Stokes-Einstein relation. 
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In the case of a simple shear along the x-axis, the expressions for the elements of 

the pressure tensor P can be obtained by substituting Eqs. and into the long 
time limit of Eq. (fT^ . For the diagonal components we obtain 



yy 



and 



pfc 



m 

/ 

m 

m 



-P 



A 



XX ^xy 



7 -Pxx[^. + l]K (26) 

(27) 

-p. (28) 
ponent different from zero is 

xy 2m ^ r-xx r-xy-xy j I 7, (29) 

where we have used Eqs. and (j77|) and defined Gq = jS^x — PxxPxyPyx — Pxxf^xyl- 
By substituting now Eqs. ((211) and into Eqs. we obtain the 

propagator 1 — M^xCi^ for the elements of the pressure tensor and from here, we may 

conclude that non-Newtonian effects appear at lower shear rates when increasing the 

radius of the particle. 

We will now proceed to analyze Eqs. (j26|) -(|29| ) by assuming 7 < In this limit, 

Eq. ()29p can be approximated up to the second order in 7 to obtain 



k 
xy 



Pxx l^x 



^xx yxy^ xy^ 



The only non-diagonal component different from zero is 



xy — 



- [riB + m] 7, 



where 



and 



1 — MxxdX-^ 



xy-^yx 



(30) 



(31) 



(32) 



In Figure 1 we present both viscosities normalized to their value at zero 7, as a function 
of the shear rate for different values of a. At low shear rates, Brownian visicosity 
(dashed-dotted and solid lines) exhibits shear thinning whereas shear thickening in the 
hydrodynamic viscosity (dotted and dahsed lines) arises at higher shear rates [21]. 

The rheological equations of state follow by substituting Eq. pO|) with (j3T| and 
(jn21) into Eqs. (j2EI) and (j2I|) obtaining 



and 



Tyk ^bT 

Kx - — P + 
m 



yy r' ^ 



yx 



M, 



xy 



''PI' 



m 



'kBT 1 
m(jo 6 , 



M^yau' 



(33) 



(34) 



from which one concludes that the lowest order contribution to the normal stress 

3 

difference is proportional to 72. Our results containing corrections to the viscosity 
have the same shear rate dependence as those obtained in Ref. [1] for simple fluids. 



Mesoscopic thermodynamics of stationary states 



8 




Figure 1. The Brownian (solid and dashed-dotted lines) and hydrodynamic (dashed 
and dotted lines) viscosities as a functions of the shear rate for different values of the 
radius of the Brownian particle. This figure shows the shear thinning and thickening 
effects. The values of the parameters we have used are /3o — 10^s~^, Dq ~ lO'^cm^s"^ 
and p ~ 0.2cTO^. 



5. Nonequilibrium equations of state for a non-ideal gas 



The thermodynamics of the gas for the case of interacting Brownian particles can also 
be formulated as in the dilute case. Within a mean field approximation in which one 
assumes that the dynamics is that of a Brownian particle moving through an effective 
medium consisting of the fluid and the remaining particles, the Fokker-Planck equation 
is 



d_ 

du 



[uP) 



knT 



m 



In/ 



(3.{u-vo)P+^i 



dP 
du 

kBT^_ dP 
du 



m 



d 






C-FP 


du 



(35) 



Interactions have been taken into account through the activity coefficient / = 
exp[^j^p~^p*"*] measuring the deviation of the Brownian gas with respect to the ideal 
case, with p*"* being the excess of osmotic pressure |18j . 

Proceeding along the lines indicated in Sec. 3, one may use Eq. in order to 
derive the evolution equation of the momentum 



P 



dv 
dt 



-p{v-vo)-(3-\/ ■F + 



knT 



m 



In/ 



Vp + pC-F, 



(36) 



where we have defined the total pressure tensor by P = P'^ The expression of 

the kinetic part of the pressure tensor is given by Pxy'^ — ~ [^b(^) + Vh{^)]i, 

with dependences of the viscosities different from those in Eq In deriving the 
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corresponding expression of P, we have assumed that the relaxation of the kinetic part 
of the pressure tensor is much faster than that of normal (structural) stresses. 

5.1. Viscoelastic effects 

As in the previous section, the transport coefficients can be identified by using the Faxen 
theorem. At higher volume fractions (p of the Brownian gas, one may assume that these 
coefficients can be expressed in terms of modified penetration lengths. 

In the creeping fiow case, the mobility coefficient can be expanded up to the first 
order in terms of the modified inverse penetration length 

a. = -r^ (no)'- , (37) 

where the scaling factor and the exponent 6^^ can, in general, be functions of the 
volume fraction 0, JH]- In E<1- (jHTj) . we have introduced the characteristic diffusion time 
r = By using Eq. (|?7j). the corresponding relation for (3^^ and e^^ can be expressed 

by m 

1 TTl 

= Po {'^ + aauj) and = -——o^Pq (1 + 2aa^) . (38) 

O KbI 

Consequently, the frequency-dependent Brownian viscosity coefficient is now given by 



Vb[^ = —rP 



l-AzHr^)'^+^Z'{roof^ , (39) 



and the hydrodynamic viscosity coefficient by 

VHico) = ^/3oaV [l + {ruY- - Aj^ {ru^f-] . (40) 
The first normal stress difference is accordingly 

P,,-FyyC:^p-'[riB{u) + rjH{uj)]j\ (41) 

and then beomes a function of frequency and volume fraction. At sufficiently low shear 
rates, the quadratic dependence on the shear rate in the first normal stress difference has 
been reported in the literature • Notice that in the ideal case, similar dependence 
can be obtained in this regime. 

5.2. Non-Newtonian effects 

In the stationary flow case, we will assume that the mobility can be expanded in terms 
of the modified inverse penetration length A^ 

A, = -r^ (tt)'^ , (42) 
A^a 

where A^ and 6^ are functions of the volume fraction (j). The mobility can be expressed 



as 



p-^ = Po\^-MaA^), (43) 
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Figure 2. Comparison of the normalized Brownian viscosity in the ideal (dash- 
dotted lines) and non-ideal cases (dashed and solid lines) as a function of the shear 
rate 7 for two values of a. In the non-ideal case, we have chosen — 0.5 and = ^ 
which correspond to a density of particles of 0.46, |18| . 



and the explicit expression of can be found by using Eq. (j2S)). In this case, the 
contributions to the viscosity are given by 

mfJo 

and 



VBh) = 1 - M,,A-' [tjY-^ + M,yMy,A-' [r^f'-'] , (44) 



Vnh) = l^oa'p [1 + M^^A-' (T7)'- + [m^ - M^yMy.,) A'^ {rj f''] (45) 
The shear rate dependent first normal stress difference is now 

P.. - Pyy ^ -^/3oaVr-V^^7^+'^. (46) 

These results show that the form of the equations of state is affected by the interactions 
by modifying the magniutde of the Brownian viscosity and the dependence of the 
exponents on the volume fraction, which is responsible for the loss of universality in 
the power law behavior of the viscosity in concentrated systems |26j . 

Figure 2 illustrates the differences in shear thinning behaviour in both the ideal 
and non-ideal cases. Since in the general case the inverse penetration length depends 
on the sum of the frequency and the shear rate |2ZI, then we have calculated the value 
of the exponent 6^ by using the results of Ref. p[B] for a density of particles of p ~ 0.46. 
From this figure and Eqs. ()44j) and ()45j) . one concludes that shear thinning effects are 
more pronounced in the non-ideal case at lower values of 7. 

6. Discussion 



By means of a reformulation of the local equilibrium hypothesis in the domain of the 
mesoscale, we have shown that a nonequilibrium thermodynamic treatment including 



Mesoscopic thermodynamics of stationary states 



11 



the presence of fluctuations and of external forces driving the system to nonequilibrium 
steady states is possible. Using this scheme, we have obtained the Fokker-Planck 
equation which in the case of an inhomogeneous bath, is coupled to the evolution 
equations of the fields characterizing the state of the bath. The evolution equations of 
the hydrodynamic fields are derived from the Fokker-Planck equation; their stationary 
solutions lead to nonequilibrium equations of state. 

We have studied, in particular, the case of a gas of Brownian particles in contact 
with a heat bath when the whole system is in motion. The ideal and non-ideal cases 
have been analyzed. In both cases, we have obtained nonequilibrium equations of state 
showing a dependence of the pressure tensor of the gas on the frequency and on the 
shear rate, thus revealing the presence of viscoelastic and non-Newtonian effects. 

In the ideal case, the non-analytical dependence on the shear rate of the corrections 
to the first normal stress difference and the viscosities coincide with those reported in 
Ref. [T] for a simple fluid. In the non-ideal case, nonequilibrium equations of state have 
been derived by assuming that the inverse penetration length is similar to that of the 
ideal case with an exponent and a scaling factor depending on the volume fraction of 
the 'gas' of Brownian particles. Our results show that, at sufficiently low shear flows, 
the first normal stress difference is a quadratic function of the shear rate thus coinciding 
with the simulations reported in Ref. [Hj. At higher shear flows, non- Newtonian effects 
are manifested through a power law dependence on the shear rate of the flrst normal 
stress difference and on the viscosities ^T] • 

Our general framework can be useful in the study of the dynamics of small systems 
subjected to driving forces or gradients for which the presence of fluctuations may play 
a signiflcant role. 
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